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detailed  and  rigorous  analysis  of  thermoelectric  transport  and  optimization  of  ZT  are  discussed. 

PACS  numbers: 


I.  POSITION-DEPENDENT  BAND  CONSIDERATION^ 

We  consider  a  sample  with  a  nonuniform  band  structure  which  is  determined  by  the  crystal  potential,  subject  to  an 
additional  internal  or  external  fields  of  force  such  as  an  electric  field,  stress,  or  temperature  gradients.  Let  (r,  k)  be 
the  band  structure  of  the  unperturbed  crystal,  where  p  is  the  band  index.  S^{ry  k)  =  S^^if^k)  -h  is  the  band 

structure  in  the  presence  of  additional  fields  of  force,  and  =  €^{f,  km)  is  the  band  edge,  where  km  are  the  wave 
vectores  corresponding  to  the  minima  or  maxima  of  the  band,  is  the  potential  energy  for  these  forces. 

Thus  the  total  energy  of  an  electron  for  band  p  is  Ep{f,  k)  ==  £73 (r,  -I-  ^(f).  We  can  introduce 

to  be  the  energy  contours  of  the  total  energy. 

For  conduction  band,  we  have, 

Ec{r,k)^Sc{r,k)^<fl{i^ 

Ec{f)^Ec{f)^n^ 

and  for  valence  band, 

E,{r,k)  =  e,{T,k)  +  n^ 

II.  WANNIER  EQUIVALENT  HAMILTONIAN^ »» 

For  a  single  band,  for  example,  the  conduction  band,  the  total  Hamiltonian  operator  can  be  replaced  by  the 
equivalent  Hamiltonian  operator. 


20050201  046 


(1) 
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where  £cir,k)  is  the  band  energy. 

The  classical  equivalent  Hamiltonian  is  introduced  by  the  correspondence  >  k, 

Hc{r,k)=Se{r,k)-\-n^^Eo{f,k) 

and  which  gives  us  the  Hamiltonian  equation,  with  momentum  p  —  hk^  as 


(2) 

(3) 

(4) 


Fn  is  the  total  force  act  on  an  electron  in  the  conduction  band. 

The  same  is  true  for  holes  in  valence  band,  by  noticing  that  the  force  act  on  holes  has  opposite  signs  as  on  electrons, 


Vp(r,k)  =  ^  =  lVi;EAr,k) 


^  ^dk 


dH„ 

dr 


=  VE,{T,k) 


(5) 

(6) 


III.  BOLTZMANN  TRANSPORT  EQUATION® 

Let  fn{f,k,t)  to  be  distribution  function  of  electrons  in  the  conduction  band.  The  Boltzman  IVansport  Equation 
(BTE)  states, 


(7[/n(r,ife,t)]  is  the  collision  integral.  Under  relaxation  time  approximation,  we  let 

r  .O  fn(r,k,i)  -  f^{r,k,t) 

- - 

where  is  the  local  equilibrium  distribution.  Thus  for  the  stationary  situation,  we  have, 

-*  /-*  7*\  «=7.f/-*r\.  ^  r  tl\  _ 

Vnir,  k)  •  V/„(r,  fc)  -I-  —  •  Vj^fnir,  k)  = - — - 

"  rn(r,  k) 

Writing  fn{r,k)  =  /°(f,  k)  +  /^(r,^),  we  find  the  perturbation  solution  for  /,J(r,fc)  as 
fn(.r,k)  = -Tnif,k) 


Vn{r,k)  ■  ^/“(r.fc)  -  .  Vi;f°{r,k) 


For  holes  in  valence  band,  the  Boltzman  TVansport  Equation  (BTE)  has  the  form, 

Pv 


.  VUiW)  +  ^  •  [- k,t)] 


Writing  /p(f,  it)  =  tf  (f,  k)  +  /i(f,  it), 


/p(r,fc)  =  --rp(f,fc) 


vAr,k)  •  V/p®(f,fc)  -  .  Vj/p"(r,£) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 
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IV.  LOCAL  EQUILIBRIUM  DISTRIBUTION,  ELECTRIC  AND  HEAT  CURRENT 


The  local  equilibrium  distribution  for  electrons  is  given  as  the  local  Fermi-Dirac  distribution, 


A  ^  r  1  ^ 

f^(r,k)=  e  + 1 


where  Cn(f)  is  the  Fermi  level,  or  the  electrochemical  potential,  for  electrons. 

The  electric  current  and  heat  current  induced  by  electrons  can  be  found  as,  respectively, 

•^n  =  -  A  / 


fn  =  --^J  d^kVr.(r,  k)f^(f‘,  k) 

/nQ  =  ^  /  d^kVr.ir,k)  [E,[?,k)  -  Cn(f)]  fl(r,k) 


“  47r3 


J  dPkT„{f,k)v„ 


{r,k)  \  vn(r,k)  ■Vf°(f,k) 


Of^  M  _  ^Ec{f;k) 


■y-X{r,k) 


dnQ  =  —^  j d^kr„{r, k)  k)  -  Cn(f)j  k)  v„{f,  k)  ■  Pf^{r, k)  - 


VEcif,k) 


(17) 


FVom  Eq.  (13)  we  have, 


V/^(r,  ^)  =  VE,(f,  -  '^C„(r-)  -  (^.(r,  k)  - 
k)  =  VgE.(f,  k)  =  §^fiv„(r,  k) 


^  “4^  / d^k^T„{,f,k)v„{P,k)vnir,k)  ■  VU^  +  [Ec(f,k)  -  Cn(r))  ^  (20) 


Rewrite  it  as. 


/n  =  an-V^-Li2.^ 

J„p  =  L«.V^-L”-^ 


=^"  =4^ 


J d^^^‘^n{r,k)v„(f,k)vn(r,k)  (2^ 

(2' 

LS'  =  ^  y  d"fc^r„(r, fc)t;„(r,  ^)^;„(f,  fc)  [Ee(r-,  fc)  -  Cn(r-)]  (2f 

^  “4^  /  ' 

is  the  electric  conductivity  for  electrons  in  the  conduction  band.  We  can  also  find  the  thermal  conductivity  as, 

JnQ  =  -«nVT  (21 

Jn=0 

^  ^  [T  22  t21  _  -1  t12-| 

^xi  —  rp  [h<n  ‘  L„  J 


(29) 
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We  can  also  define, 


=  S„VT 

®  /n=o 

s 

On  —  jt 


For  holes, 


/p(^*:)  =  l-/n(^>*)  =  l-  e  +1  =  ' 


e  ““  +1 


where  Cp(^  is  the  Fermi  level,  or  the  electrochemical  potential,  for  holes. 

The  electric  current  and  heat  current  induced  by  electrons  can  be  found  as,  respectively, 

4  ^  y  d^kvpir,  k)fp (r,  k) 

JpQ  =  -^  j d^kvp{f,k)  [E„(f,k)  -  <p(f)j  f^(f,k) 


j^pQ  =  -^  I  k)  [£?„(r-,  fc)  -  CpCO]  Vp(f,  k)  Vp{f,  k)  ■  k)  -  k) 


FVom  Eq.  (32)  we  have. 


Vf°{f,k)  =  VE,{f,k)  -  V<p(f)  -  [E,(r,k)  -  Cp(f)) 


~  4^  / fc)  -  Cp(f))  ^ 

~  4^  /  ~  ■  ^Cp(»0  +  (B„(f,  fe)  -  Cp(iO)  ^ 


Rewrite  it  as, 


/pg=L?^.V%-L”.$: 


Lp''  =  “Lp^  =  J  d^k-^^Tp(r,k)Vp{r ,k)vp{f ,k)  ^£„(r,fc)  —  Cp(0] 
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cTp  is  the  electric  conductivity  for  holes  in  the  valence  band.  We  can  also  find  the  thermal  conductivity  as, 


— /CpVT 


We  can  also  define, 


SpVr 


S  -a 

Op  —  CTp  ^ 

For  lattice  heat  flux  density,  we  assume  that 

Jlq  =  — V  T , 

where  Ki,  is  the  lattice  thermal  conductivity  assumed  to  be  known  from  other  experimental  considerations. 


V.  SPHERICAL  BAND,  LONGITUDINAL  ACOUSTIC  WAVE 

Assume  we  have  spherical  parabolic,  effective  mass  band  structures  centered  at  k  —  (0,0,0). 
For  electrons  in  conduction  band, 

m*  (r)  is  the  coordinate  dependent  effective  mass.  Thus  we  have, 

Mr,  it)  = 

m*(r) 

We  further*’ assume^  the  longitudinal  acoustic  dispersion  relations  for  phonons,  =  voq,  so  that, 

A  TnOkif)  T„o(f) 


Tn{r,k)  =  T„(f,k)  = 


where  r„o*(f)  = 

Component  of  the  electric  conductivity  has  the  form. 


/fe(f',fe)-£c(f) 


jSkj£M?,k)  j 


All  the  components  in  the  above  integral,  except  for  ka  and  k^,  are  spherically  symmetric,  thus  (t„  is  a  scalar  instead 
of  a  tensor  and  we  have, 


i,<rn)a0  =  0,  (a^p) 

On  =  («^n)n  =  {On)22  =  (<rn)33  = 


(gn)ll  +  (<rn)22  +  («^n)33 


f,3rdf^ 


For  the  same  reason,  L^^,  and  Lj^  are  all  scalars, 


Lr  =  -t 


(60) 
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e^{r,k)  -  £c{f)  _  Ec{f,k)  -  Ec{f) 


Vnif^  = 


kgT 
Ec(f)  - 
ksT 


M  =  J_M 

dEc  ksT 


Let  J  dk^  —  J  dkk^  Jq  d9  sin  0  Jq^  d<j>  and  integrate  over  0  and  <j>,  noticing  that  k  =  yj  -^Vk^y/E^  and  Eq.  (54), 

(t)*  r^(m,T)f4e,e„^£  (68 


Ll^  =  Li?  = 


^  (^)  *  rUmsT?  j"  d£„£^  {£„ 


Consider  integral,  where  9{Sn)  is  an  arbitrary  function  with  g{0)  =  0, 


Then  we  have 


fir  n(p  —  \f^{F  \n(F  d£  d£  fO 

[/„(£„)fl(fn)]|^  d£„  fn-  d£„  /„ 

(¥)* f<‘S,ae..n 

3 

L]?  =  (^)  "  rM^ikeT)^  j^  d£„f°{£„,f)  [2f„  - 

Lf  =  Li^  =  rUmsT)^  d£„f°{£„,f)  [2£„  -  7^1 

Li?  =  (^)  '  r„o{f^ikBTf  jy  d£nf^{S.,i^  [3£?  -  4»j„5„  +  vl] 


Introduce  Fermi-Dirac  function  Fi(T/),  which  is  defined  as  Fi{r])  =  Jq 


D  e?^-»T)+l 


.  Thus, 


Stt^  m: 


3 

(^)  '  [2F^(Th)  -  VnFoiVn)] 

3 

=  (^)  ^  [2Fi(7,„)  -  VnFo(rin)] 


r  22  _  _ f 

”  m*  \  h 


)(f)(A:BT)  [3^2(7^71)  ^VnFl{T}n)  ~i~ ’Hn^oiVn)] 
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and  then, 


5n(r)  = 


ks  2Fi  {T]n)  -ynPoirtn 


ffnT  e  Foirin) 

<^nT 

Kn{f)  X,22_i^  Ub/  7 


=  e  V  1  ^o^(^n) 

Ub  /  T  3Fo(»j„)F2(»J„)  -  4F2(,,„) 


For  holes  in  valence  band, 


mp(f)  is  absolute  value  of  the  coordinate  dependent  effective  mass.  Thus  we  have, 

We  again  assume  the  longitudinal  acoustic  dispersion  relations  for  phonons,  cjqi  —  voq,  so  that. 


Tp(f,k)  =  Tf,{f,k) 


'rpOfc(f) 


l€p{f)-Sp{r,k) 


where  r„op(r)  =  and  rpo(f)  = -^TpOfcW. 

For  the  same  reason  as  for  electrons  in  conduction  band,  cTp,  Lp^,  Lp^  and  Lp^  are  all  scalars. 


127r^ 


-C,(f)] 


7;p(f)  = 


£„{f)  -  Sy{r,k)  _  Eyjf)  -  Ev{f,k) 
ksT  ksT 

K(r-)-Cp(0 


then  we  have. 


_ i_^ 

ksT  dSp 


(78) 

(79) 

(80) 

(81) 

(82) 

(83) 

(84) 

(85) 

(86) 

(87) 

(88) 

(89) 


(90) 
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Let  /  dk^  —  f  dkk^  Jq  d9  sin  0  Jq^  d(l)  and  integrate  over  6  and  0,  noticing  that  k  =  y^^y/ksT^/E^  and  Eq.  (82), 


Stt^  m*  \  h 


=  ‘^p^p(»7p-4)^ 


r21  _  rl2  _  ^ 

-  -t-n  -  3^2^; 


1  1  /2m!\^ 


P  “  35r2  1 


9m*  \  5  /•°°  of** 

-/]  rAmBT?  j^  ^P^P(»?P-^P)^^ 


Consider  integral,  where  ^(fp)  is  an  arbitrary  function  with  p(0)  =  0, 


f°°  Af  n(F  l^'^P  —  \f^(F  \n(£  'll  —  f°°  d£  —  _  f°°  A£ 

Jo  -  [A.(^p)ff(^p)J  ^  ds^  u-  ^p  QS^  f, 


^si^p)  fO 


Thus, 

-  1  /im*  \  i  /“” 

=  d^p(»fc-2fp)/«(£:p,rl  (97) 

'3 

L2i  =  -Li2  =  ^i^^)%^(f)(A:BT)2^“dfp(,,p-25p)/«{£,,r)  (98) 

if  =  -3^;^  (^)  ’  Tpo(r-)(fcBT)2 d£,  [Zel  -  4np£p  +  V^]  f°{£p,^  (99) 

Thus, 

3 

<Tp  =  ^;^(^)%po(»0(feBT)f'o(7;p)  (100) 

if  =  (^)  ’  r^ornksTf  [2Fi(7,p)  -  7,pFo(r,)]  (101) 

♦  ^ 

if  =  -Lf  =  - 3^ (^)  ’  r^(^0(fcBT)=*  [2Fi(,,,)  -  r^Fo(v)]  (102) 

3 

if  =  -^;^(^)’rpo(r-)(*Br)^[3F2(»?p)-%Fi(,^)  +  r,^Fo(7,p)]  (103) 


and  then, 

q  f'^—  (^p)  Vp^oi'np) 

apT  e  Fo(77p) 

crp(f)  (TpT  _  /"  e  y  1  F^(r?p) 

«p(r-)  L22_ii!£li  Ub/  T3Fo(»?p)F2(7/p)-4F2(,ft,) 

P  tTp 


(104) 

(105) 
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VI.  NONDEGENERATE  BOLTZMANN  APPROXIMATION 


In  the  nondegenerate  approximation,  Ec{f,  >  Ec{^-C,n{^  >  and  (^p{f)-Ev{r,  k)  >  Cpi^^^vi‘0  ^ 

A:flT,the  Fermi-Dirac  distribution  function  goes  to  classical  Boltzman  distribution, 


E..(f’.E)-<p(*0 

fpif,k)-*e  '‘bt 

(106) 

and 

Jo 

(107) 

That  is, 

roo 

Foiv)->  /  =  e” 

Jo 

poo 

Fiiv)-*  /  £e-^^-^U€  =  e^ 

Jo 

poo 

F2(»?)-»  /  £^e-^^-^'>d£  =  2e’> 

Jo 

(108) 

Thus, 

(109) 

(110) 

-  -3i;4  (^) '  -  ’>•1 

(111) 

i?  =  3^5  (t)  ' 

(112) 

(113) 

(114) 

=  -3^;^  (^)^p.(rl(^.r)V.  p  -  »| 

(115) 

if  -  -5^^  (^)  '  VimiTfe"-  [6  -  4%,  +  vf 

(116) 

and  then. 

„  2kB(.  Vr.\  2kB  f 

e  V  Y)-  e  \}^  2kBT{f)  ) 

(117) 

(r„(f)  a„r  V  ^ 

K„(f)  i22  2Ub/  T 

n  (Tn 

(118) 

(119) 

<^pT  __ 

Kp{r)  iji2  2  \A:b/  T 

P  <Tp 

(120) 
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VII.  TWO  BAND  -  CONDUCTION  BAND  AND  VALENCE  BAND  -  MODEL 
Letting  it  follows  in  a  straightforward  fashion  that  the  two-band  dynamical  variables  become 

J  =  /„  -b  jp  =  [On  +  CTp)^  -  ^ 

where  cr  =  cTn  +  cTp,  and  L02)  =  £,0^1.  From  this,  we  obtain 

g-n  5n+0‘p  «Sp 


(ffn+ffplT  (Tn  -l-CTp 


/o  =  /ncj  +  /pg  =  +  Lf’}  ^  ^ 

where  and  -h  .  Also,  we  acquire 


»  -  j,  (iM -  fy)  -  j,  ^ig“' + 4“>  - 
.  I  (tk  +  +  Tk,  +  a-jg  -  'I 

T  <T„  Po^  »-„+<T,  y 

=  fc„  -b  fcp  -b  ^  {T5„z,<i*^^  -b  rSpLj?^’  -  rs(42^>  +  42^>)} . 


K  =  kn  +  kp  +  {S„^  S)  Ll^i)  -b  (5p  -  5)  Lf  . 


VIII.  SEEBECK  COEFFICIENT;  FIGURE  OF  MERIT^ 

For  a  given  carrier  band,  the  current  is  given  by  Jn  in  eq.  (22)  and  Jp  in  eq.  (41),  respectively.  Then,  since  the 
Seebeck  coefficient  is  determined  from  the  open  circuit  condition  that  Jn  or  Jp  =  0,  it  follows  that 

(IM) 

where  m  represents  n  or  p  from  eqs.  (22),  (41),  respectively.  Here,  Cm(^  is  the  electro-chemical  potential,  given  by 

Cm{f)  =  Ef„  +q<f>{r),  (127) 

where  <t>{f)  is  the  external  potential  across  the  sample  and  q  is  the  charge  of  the  carrier.  Therefore,  after  integrating 
across  the  sample  of  length  L,  one  obtains 

^(L)-.^(0)  =  {£lF„(0)-£F„(l^)}  +  /o''5„Vrdx,  (128) 

where  Sm  —  T)  is  defined  in  eqs.  (31), (50)  for  m  =  n,p,  respectively.  For  aslab  in  which  £?Fm(0)  =  Ef^(L), 

we  see  that  eq.  (128)  simplifies  so  that  in  one  dimensional  case,  we  get 

=  <I,{L)  -  m  =  qf^^Sm^dx.  (129) 

Therefore,  from  eq,  (129)  we  can  see  that  the  voltage  developed  across  the  slab  depends  upon  the  quantities  5m  (x)  = 
Lm^\x)/{ayn{3:)T{x))  and  and  am  have  already  been  calculated  for  conduction  and  valence  band  cases, 

and  for  degenerate  and  non-degenerate  Boltzmann  approximations.  As  for  T(x),  given  that  the  temperature  at  the 
ends  of  the  slab  are  fixed  at  hot  and  cold  values  to  provide  a  temperature  gradient  across  the  sample,  the  value  of 
the  temperature  within  the  sample  will  vary  with  position  since  T{x)  will  classically  depend  upon  the  solution  of  a 
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boundary  value  problem  for  thermal  diffusion.  The  classical  solution  has  been  well  documented,  and  in  fact,  in  the 
open  circuit  condition,  it  is  well  known  that  T{x)  varies  linearly  across  the  slab  as 

Tix)  =  TiQ)  +  ^x,  (130) 

where  T(0)  is  the  cold  temperature  and  T(L)  is  the  hot  temperature,  and  AT  =  T(L)  -r(0),  Using  T{x)  of  eq.  (130) 
in  eq.  (129)  therefore  gives  the  simple  result  that 

(131) 


That  is,  the  Seebeck  coefficient  is  the  spatial  average  of  Sm{^)  over  the  slab. 

We  have  evaluated  Sm  =  Lm^\x)/{am{x)T{x))  for  the  spherical  band  model  with  spatially  varying  band  edge  and 
longitudinal  acoustic  wave  phonon  approximation  in  Section  V.  Sm  has  been  evaluated  as  a  function  of  position  for 
both  degenerate  Fermi-Dirac  [eqs.  (78),  (104)  for  n,p,  respectively]  and  non-degenerate  Boltzmann  thermodynamic 
limits  (eqs.  (117),  (119)  for  n,p,  respectively).  For  simplicity  of  discussion,  we  present  the  result  for  the  non-degenerate 
case  for  electrons  so  that  eq.  (131)  becomes 

^  =  -r  So"'  ^  0  +  ^  (132) 

or 

s  =  (2  +  ifo^  dx)-  (133) 

Thus,  the  Seebeck  coefficient  is  enhanced  by  the  addition  of  the  spatial  average  of  spatially  dependent  Boltzmann 
energy  ratio  {Ec  —  Cn)/i^BT)  directly  in  the  Boltzmann  limit.  If  we  were  to  evaluate  Sn  for  electrons  in  the  Fermi- 
Dirac  limit,  we  would  have  to  use  in  eq.  (78)  which  would  then  require  an  integral  over  Fermi-Dirac  integral 
functionals  of  {Ec  -  Cn)/{^BT)  in  a  much  more  non-linear  dependence  on  the  spatially  varying  Boltzmann  energy 
ratio. 

For  simplicity  in  discussing  the  estimate  of  the  enhancement  factor,  we  assume  that  Cn(a;)  =  Cn(0)  a  constant.  We 
note  that  if  the  Boltzmann  factor  were  a  simple  constant  in  total,  the  enhancement  would  simply  result  in  an  increase 
in  the  Seebeck  coefficient  from  the  free  electron  value  to  the  semiconductor  value.  In  fact,  in  enhancing  the  Seebeck 
coefficient,  we  simply  want  to  maximize  the  area  under  the  curve  of  the  Boltzmann  energy  ratio  in  eq.  (133).  However, 
this  dependence  upon  the  Seebeck  coefficient  needs  to  be  balanced  with  the  other  transport  parameters  that  enter 
into  the  well-known  figure  merit,  Zm{x)  which  is  given,  for  a  given  carrier  type,  as 


Z  r(22)  "(12)  ^(21) 

‘■m  +kL 


(134) 


Here,  there  dependence  of  .  and  upon  the  spatially  dependent  Boltzmann  factor  has  been 

calculated  for  degenerate  and  non-degenerate  statistics  in  Section  V,  VI,  respectively;  also  is  the  lattice  thermal 
conductivity.  For  the  case  of  non-degenerate  statistical  analysis,  we  see  from  eqs.  (109)-(112)  for  electrons  that 
=  L^m^  Can  be  expressed  in  terms  of  so  that  Zn  in  eq.  (134)  can  be  greatly  simplified  as 


ZTn{^)T{x)  2  a^{x)  f{x)  +  /cl  2  T(x))  ’ 


(135) 


where,  from  eq.  (109),  (t„  is  given  by 

(^)'^"^"o(x)(fcBT(x))e-(^«-<«)/(‘‘^W).  (136) 

The  dependence  of  all  the  thermodynamic  variables  on  the  spatially  dependent  Boltzmann  factor,  even  though  in¬ 
troduced  explicitly  here  through  the  approximate  form  of  the  relaxation  approximation,  is  resonant  with  Landauer’s 
original  concept  of  “motion  out  of  noisy  states”  ®  where  non-uniform  temperature  effects  drive  the  system  dynamics 
-  in  this  analysis,  it  is  a  non-uniform  Boltzmann  factor  which  drives  the  system,  where  the  energy  band  edge  is 
tailored  and  fixed,  and  the  non-uniform  temperature,  apart  from  the  fixed  hot  and  cold  ends,  is  transferring  energy 
dynamically.  Therefore,  even  though  {Ec  -  Co)/{kbT{x))  can  be  chosen  to  enhance  5„(x)  over  the  slab  of  material, 
the  spatial  dependence  of  the  figure  of  merit  will  depend  upon  the  term  kL/{(Xn  T{x))  in  eq.  (135)  in  a  very  sensitive 
way;  in  fact,  the  model  dependence  of  (T„  in  eq.  (136)  shows  that  the  conductivity  will  depend  exponentially  on  the 
spatial  Boltzmann  factor. 
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As  a  few  illustrative  examples  of  Ec{x)  to  estimate  S  in  eq.  (132),  we  consider  the  following  simple  cases,  always 
assuming  that  T{x)  =  T(0)  +  ^  x,  with  AT  =  T{L)  -  T(0),  and  C  =  C(0)  =  Co- 
Case  1:  (2:)  =  £'c(0),  a  constant  band  edge. 


=  i  lo  ^  (137) 

Case  2:  Ec{x)  =  Ec(0)  +  ^  x,  a  linear  spatial  band  edge. 

=  i  /cf  dx  ^  In^  +  ^  (1  -  ^  1„|^)  ,  (138) 

Case  3:  Ec{x)  =  Ec{0)-^Yln^o  8{x  —  na),  aperiodic  but  asymmetrical  spatial  sawtooth  band  edge;  here,  S(x-na) 
defines  a  unit  cell  (see  sketch)  of  the  periodic  band  edge,  where 


£{x  —  n  a)  = 


^{^-(n-1)}  for  (n  -  l)a  <  X  <  (n  —  l)a  4*1 

^  (_^  4.  li)  for  (n  -  l)a  H-  ra  <  X  <  na 
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Comparative  inspection  of  eqs.  (142),  (143)  indicates  that  A/3(l)  <  A/3(l/2)  so  that  for  all  other  equal  parameters, 
the  symmetric  periodic  band  edge  produces  a  larger  area  under  the  curve  of  eq.  (140)  than  the  asymmetric  band  edge; 
thus  the  enhancement  to  the  spatially  averaged  open-circuit  Seebeck  coefficient  would  be  larger  for  the  symmetrical 
(r  =  1/2)  band  edge  case.  However,  it  should  be  noted  that  even  though  a  given  choice  of  spatial  energy  band 
configuration  may  give  rise  to  a  relatively  large  Seebeck  coefficient,  the  more  significant  consideration  is  to  choose  an 
energy  band  configuration  which  optimizes  the  figure  of  merit  ZT.  In  this  case,  the  physics  is  not  determined  by  open 
circuit  conditions,  as  is  done  in  evaluating  the  Seebeck  coefficient,  but  is  determined  generally  under  non-equilibrium, 
current  driven  conditions.  Thus,  an  effort  should  be  pursued  to  examine  the  role  of  band  engineering  on  the  relevant 
thermoelectric  transport  properties,  and  to  properly  address  the  ability  to  enhance  and  optimize  the  thermoelectric 
figure  of  merit  based  on  a  basic  variational  principle  for  ZT. 


IX.  VARIATIONAL  PRINCIPLE  FOR  ZTy  THE  FIGURE  OF  MERIT 


As  noted  in  eqs.  (134),  (135),  the  well-known  figure  of  merit  for  a  given  carrier  type  is  generally  given  by 


ZmT  = 


+  • 


(144) 


In  a  very  nice  variational  analysis  developed  by  Nishio  and  Hirano,®  they  treated  Sm  and  as  functionals  of 
differential  conductivity,  t7£)(6),  where  aoie)  is  known  and  is  defined  as  the  integrand  of  total  conductivity 


cTm  =  /<r£)(e)de.  (145) 

The  transport  variables  and  Km  can  be  expressed  in  terms  of  <tx)(c)  as 

-  C)<^£>(e)  (146) 

and 

K„.  =  0^<^D(e)  -  Ta^Sl ,  (147) 

where  C  is  the  chemical  potential,  and  Kl  the  lattice  thermal  conductivity,  is  considered  constant.  Thus,  ZmT  in 

eq.  (144)  can  be  viewed  as  a  functional  of  a d  alone,  provided  ATi,  is  a  constant.  Thus,  a  small  change  in  (Td{^) 

functionally  gives  rise  to  a  change  in  Z  through  the  functional  derivative  expression 

6Z„,  =  Zl<Jm  +  Sa„]  -  Z[a^]  =  fde^^5(D(e).  (148) 


Nishio  and  Hirano®  have  shown  that  this  functional  derivative  in  eq.  (148)  yields  a  maximum  for  eq.  (144)  in  an 
energy  range  given  by 

^  =  -"I-  (i  +  ^)  (i±7rbr)-  (1^9) 

This  suggests  that  it  is  possible  to  improve  the  figure  of  merit  of  thermoelectrics  by  tailoring  the  transport  energetics 
by  band  engineering  the  carriers  so  that  they  transport  across  the  thermoelectric  device  in  an  energetically  favorable 
range. 

The  variational  principle  developed  by  Nishio  and  Hirano  is  very  illuminating,  but  it  only  applies  to  situations 
where  the  lattice  thermal  conductivity  is  constant.  A  more  useful  form  of  the  variational  principle  would  allow  the 
lattice  thermal  conductivity  to  vary  as  a  subsidiary  condition  through  a  functional  dependence  on  the  conductivity. 
It  is  interesting  to  note  that  the  expression  for  ZmT  in  eq.  (135)  can  be  written  in  terms  of  a  calculated  function  of 
transport  variables  which  can  all  be  expressed  in  terms  of  the  conductivity  Therefore,  ZmT 

would  be  a  function  of  {Ec  —  Co)/kT{x)y  the  spatially  dependent  relaxation  time  rno(x),  and  the  spatially  dependent 
eflFective  mass  m*(x).  In  plotting  as  a  function  of  (Ec  —  C)/^Ty  one  would  find  that  ZmT  peaks  for  a  given 
thus  indicating  that  ZmT  is  amenable  to  optimization  with  respect  to  the  Boltzmann  ratio  parameter  in  a(x). 
This  optimization  scheme  is  given  physical  meaning  based  on  the  variational  principle  discussed  herein. 


X.  DISCUSSION  FOR  FOLLOW-ON  AND  CONCLUSION 

In  this  study,  we  have  explored  the  influence  of  spatially  graded  energy  bands  on  the  thermoelectric  properties  of 
thin  film  semiconductors.  For  an  assumed  linear  spatial  temperature  variation  across  the  sample  during  open  circuit 
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conditions,  the  Seebeck  coefficient  was  determined  and  was  shown  to  be  enhanced  by  the  addition  of  a  term  which 
depended  analytically  upon  the  Boltzmann  factor  comprised  of  the  “band-engineered”  energy  band  edge  divided  by 
the  spatially  dependent  temperature  across  the  sample.  Estimates  of  the  enhanced  Seebeck  coefficient  for  a  symmetric 
versus  asymmetric  periodic  band  edge  indicate  that  the  symmetric  band  edge  gave  rise  to  a  larger  Seebeck  correction; 
however  we  noted  that  the  more  significant  consideration  for  ultimate  thermoelectric  optimization  is  the  figure  of 
merit,  which  is  determined  under  current-driven  conditions  -  here  we  feel  that  the  asymmetric,  periodic  band  edge 
will  have  a  dominant  role  in  driving  the  non-equilibrium  carrier  dynamics  due  to  its  non-linear  influence  upon  the 
conductivity.® 

Overall,  the  study  revealed  the  potential  for  optimizing  the  thermoelectric  figure  of  merit  through  band-gap  en¬ 
gineering.  A  more  detailed  and  rigorous  analysis  of  thermoelectric  transport  under  graded  band-gap  conditions  are 
required,  and  a  more  comprehensive  approach  to  the  variational  optimization  of  ZT,  including  the  integration  of 
thermal  lattice  dynamics,  is  necessary  to  explore  the  physical  frame  work  for  uncovering  the  potential  of  enhanced 
thermoelectricity  in  thin  film  semiconductors.  This  more  detailed  and  rigorous  analysis  moves  along  the  following 
road  map. 

•  formulate  the  transport  equations  with  a  careful  analysis  of  the  spatial  dependence  of  the  energy  bands;  this  means 
developing  a  Wannier-equivalent  theorem  for  spatially  varying  band  edges  from  the  Wigner-Boltzmann  picture,  and 
deriving  the  correct  equation  of  motion  for  the  Boltzmann  drift  term. 


•  develop  the  full  Boltzmann  formulation  for  carriers  in  bands  subject  to  electron-phonon  interactions,  and  electron- 
mediated  phonon  dynamics.  One  needs  the  electron  Boltzmann  equation  with  electron-phonon  scattering;  in  addition, 
one  needs  the  phonon  Boltzmann  equation  with  electron-phonon  scattering.  Both  equations  are  solved  simultane¬ 
ously  with  inhomogeneous  energy  band  edges  and  doping,  and  then  caiculates  the  Seebeck  coefficient  and  ZT  using 
appropriate  boundary  conditions. 

•  explore  the  variational  principle  for  Z  T  optimization  including  both  fixed  and  varying  lattice  thermal  conductivity 
as  a  subsidiary  condition. 


•  overall  objective  -  fully  explore  the  role  of  band-gap  engineering  in  enhancing  the  thermoelectric  properties  of 
semiconductors 
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